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1. Introduction 



The discovery of quasicrystals |jj has focused considerable interest on quasiperiodic or, 
more generally, aperiodic systems [Q]. In the field of critical phenomena, due to their 
intermediate situation between periodic and random systems, aperiodic models have 
been intensively studied (for a review, see 0]). Furthermore, aperiodic multilayers are 
experimentally feasible and should build a new class of artificial structures exhibiting 
interesting bulk and surface properties. Although aperiodic superlattices have already 
been worked out by molecular beam epitaxy @, nothing has been done experimentally 
up to now from the point of view of critical phenomena. In the perspective of possible 
future experimental studies in this context, it seems an interesting and challenging 
problem to complete our understanding through a mean field theory approach. Surface 
critical behaviour has indeed been intensively investigated on the basis of the Ginzburg- 
Landau theory || in the seventies ||. This led to a classification of the transitions 
which may occur at the surface and to the derivation of scaling laws between surface 
and bulk critical exponents (for a review, see ||). These early papers are known as 
an important stage in the further developments of surface critical phenomena. 

Seen from the side of critical phenomena, the universal behaviour of aperiodically 
perturbed systems is now well understood since Luck proposed a relevance-irrelevance 



criterion 0, ITOf. The characteristic length scale in a critical system is given by the 



correlation length and as in the Harris criterion for random systems flll |, the strength 
of the fluctuations of the couplings on this scale determines the critical behaviour. 
An aperiodic perturbation can thus be relevant, marginal or irrelevant, depending on 
the sign of a crossover exponent involving the correlation length exponent v of the 
unperturbed system and the wandering exponent uj which governs the size-dependence 



of the fluctuations of the aperiodic couplings |T2| . In the light of this criterion, the 



results obtained in early papers, mainly concentrated on the Fibonacci |T3| and the 
Thue-Morse [Q sequences, found a consistent explanation, since, resulting from the 
bounded fluctuations, a critical behaviour which resembles the periodic case was found 
for the Ising model in two dimensions. 

In the last years, much progress have been made in the understanding of the 
properties of marginal and relevant aperiodically perturbed systems. Exact results 
for the 2d layered Ising model and the quantum Ising chain have been obtained with 
irrelevant, marginal and relevant aperiodic perturbations |15, 16]. The critical behaviour 
is in agreement with Luck's criterion leading to essential singularities or first-order 
surface transition when the perturbation is relevant and power laws with continuously 
varying exponents in the marginal situation with logarithmically diverging fluctuations. 
A strongly anisotropic behaviour has been recognized in this latter situation []I7, [18[ . 
Marginal surface perturbations have also been studied with the Fredholm sequence [[0| 
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and conformal aspects have been discussed [f20|j . 

In the present paper, we continue our study of marginal sequences. The case of 
the Fibonacci sequence, which leads to irrelevant behaviour in the Ising model, should 
exhibit non universal properties within mean field approach according to the Luck 
criterion and it has not yet been studied in this context. The article is organized 
as follows: in section |2], we present the phenomenological Ginzburg-Landau theory on a 
discrete lattice with a perturbation following a Fibonacci sequence and we summarize 
the scaling arguments leading to Luck's criterion, then we discuss the definitions of 
both bulk and surface thermodynamic quantities. We consider magnetic properties in 
section || Both bulk and surface quantities are computed numerically, leading to the 
values of the corresponding critical exponents. In section f|, we discuss the thermal 
properties and eventually in section [5], we discuss the upper critical dimension of the 
model. 



2. Discrete Ginzburg-Landau equations for a Fibonacci aperiodic 
perturbation 

2.1. Landau expansion and equation of state on a one- dimensional lattice 

Let us first review briefly the essentials of the Ginzburg-Landau theory formulated 
on a discrete lattice. We consider a one-dimensional lattice of L sites with a lattice 
spacing i and free boundary conditions. The critical behaviour would be the same 
as in a d— dimensional plate of thickness Li with translational invariance along the 
d — 1 directions perpendicular to the chain and extreme axial anisotropy which forces 
the magnetic moments to keep a constant direction in the plane of the plate. We 
investigate the critical properties of an aperiodically distributed perturbation within the 
framework of a 4 phenomenological Landau theory The underlying assumption 



in this approach is based on the following expansion of the bulk free energy density 

hi*,} = + \^ ~ Hfr + \c (h±l^ 2 , (l) 

where the aperiodic perturbation of the coupling constants is determined by a two-digits 
substitution rule and enters the (ft 2 term only. A dimensional analysis indeed shows that 
the deviation from the critical temperature, /i, is the relevant scaling field which has to 
be modified by the perturbation. The free energy of the whole chain is thus given by 

%]=EMW. ( 2 ) 

3 

and the spatial distribution of order parameter satisfies the usual functional 
minimization: 

6F[<f>j] = F[(j)j + 6^} - F[<f>j] = 0. (3) 
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One then obtains the coupled discrete Ginzburg-Landau equations: 



Hih + gcf)] -H- ^(<Pj+i - 20, + = 0. (4) 
The coefficients \ij depend on the site location and are written as 

fij = k B T - (zJ - Rfj) =a Q (l-- + rf^j , (5) 

where J is the exchange coupling between neighbour sites in the homogeneous system, 
z the lattice coordination and fj the aperiodically distributed sequence of and 1. The 
pref actor oq = k B T is essentially constant in the vicinity of the critical point, and the 
temperature 9 is normalized relatively to the unperturbed system critical temperature: 
9 = k B T/zJ. In the following, we will also use the notation fi = 1 — 1/9. In order to 
obtain a dimensionless equation, let us define 4>j = rrij \Jao/g leading to the following 
non-linear equations for the m/s: 

(/i + rfj)rrij + w? j -h- (m j+ i - 2rrij + m j _ 1 ) = 0, (6) 
with boundary conditions 

(/i + rfxjmi + m\ — h — (m 2 — 2m\) = 0, (7a) 

(A* + r fL)m L +m\-h- {-2m L + m L -i) = 0. (7b) 



Here, the lengths are measured in units £ = yc/ao and h = H^g/al is a reduced 
magnetic field. 

One can point out the absence of specific surface term in the free energy density. 
The surface equations for the order parameter profile simply keep the bulk form with 
the boundary conditions m = rriL+i = and our study will only concern ordinary 
surface transitions M. 

2.2. Fibonacci perturbation and Luck's criterion 

The Fibonacci perturbation considered below may be defined as a two digits substitution 
sequence which follows from the inflation rule 

0^5(0) =01, 1^S(1)=0, (8) 

leading, by iterated application of the rule on the initial word 0, to successive words of 
increasing lengths: 



1 

010 (9) 
10 1 v ' 

10 10 10 
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It is now well known that most of the properties of such a sequence can be characterized 
by a substitution matrix whose elements My are given by the number n^' of digits of 
type i in the substitution S(j) M, O. In the case of the Fibonacci sequence, this yields 



n* (0) 4^ \ Ml 



M=[ ^(o) ^5(1) | = ( , .1- do) 



1 



The largest eigenvalue of the substitution matrix is given by the golden mean Ai = 1+ 2 V ^ 
and is related to the length of the sequence after n iterations, L n ~ A™, while the 
second eigenvalue A 2 = — 1/Ai governs the behaviour of the cumulated deviation from 
the asymptotic density of modified couplings Poo = 1 — ^ + { - 

E(/i - /) =n L - PooL H A 2 |»~ (A*)", (11) 
i=i 

where we have introduced the sum Ul = /j an d the wandering exponent 

When the scaling field p is perturbed as considered in the previous section, 

Hi = aoin + rfj) , (13) 
the cumulated deviation of the couplings from the average at a length scale L 
— 1 L 1 

Sfi(L) = - 22{pj - fi) = T a r(n L - p^L) (14) 

L 3=1 L 

behaves with a size power law: 

JJ1(L) ~ L- 1 , (15) 

and induces a shift in the critical temperature St ~ ^ w ^ 1 to be compared with the 
deviation t from the critical temperature: 

**„ r M«-i)+i). ( 16) 

This defines the crossover exponent = — 1) + 1. When = 0, the perturbation 
is marginal: it remains unchanged under a renormalization transformation, and the 
system is thus governed by a new perturbation-dependent fixed point. 

A perturbation of the parameters g or c entering the Landau expansion ([[]) would 
be irrelevant. 
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2.3. Bulk and surface thermodynamic quantities 

In the following, we discuss both bulk and surface critical exponents and scaling 
functions. We deal with the surface and boundary magnetizations m s and m 1; surface 
and boundary susceptibilities Xs and \ii an d surface specific heat C s . All these 
quantities can be expressed as derivatives of the surface free energy density f s |§ (see 
table I). 

Table 1. Bulk and surface thermodynamic quantities in terms of the bulk /t> and 
surface f s free energy densities, h and hi are bulk and surface magnetic fields 
respectively and t is the reduced temperature. 





magnetization 






susceptibility 






specific heat 




bulk 


surface 




bulk 




surface 




bulk 


surface 




m b = 


dh " ts 

TOi = — 


dfs 
dh 
dfs 
dhx 


Xb = 


d 2 h 


Xs = - 

Xi = - 
Xn = - 


d 2 f s 

Uh T 

d 2 .fs 
dhdhi 
d 2 fs 
~dhf 


c b = 


9 2 h r _ 


d 2 f s 



While there is no special attention to pay to these definitions in a homogeneous 
system, they have to be carefully rewritten in the perturbed model that we consider 
here. First of all, we shall focus on local quantities such as the boundary magnetization 
mi or the local bulk magnetization m(„_!), defined, for a chain of size L n obtained 
after n substitutions, by the order parameter at position L n _i. This definition leads to 
equivalent sites for different chain sizes (see figure [l]). 

fj 1 1 1 1 00101001 010 
o«oo*o«oo*oo*o»oo«o«o 

J 1 2 3 4 13 21 



m 1 m (n-l) m L 

Figure 1. Fibonacci chain of 21 sites. The local bulk magnetization, for a chain of L n 
sites obtained after n iterations of the substitution rule is computed on the site 
here site 13. 

In addition to these local quantities, one may also calculate both surface and 
mean bulk magnetizations (m s and m& respectively), which should be interesting 
from an experimental point of view since any experimental device would average any 
measurement over a large region compared to microscopic scale. In order to keep 
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symmetric sites with respect to the middle of the chain, and to avoid surface effects, the 
mean bulk magnetization m b is defined by averaging over L n _ 2 sites around the middle 
for a chain of size L„. 



E 



rrii 



L n-2 j£Ln _ 2 



;i7) 



We checked numerically that one recovers the same average as for a chain of size L n _ 2 
with periodic boundary conditions. Following Binder || , for a film of size L n with two 
free surfaces, the surface magnetization is then defined by the deviation of the average 
magnetization (rrij) over the whole chain from the bulk mean value: 



ni 



- - 2 I m b 




A graphical description can be found in figure 0. 




Figure 2. Typical shape of the order parameter profile for a perturbed system, 
showing the boundary and local bulk magnetizations mi and mi n -V)i and the average 
values m& and (rrij) 

In the following, we shall use brackets for the averages over the finite system, taking 
thus surface effects into account. In the same way, the bulk free energy density in table [I] 
has to be understood as: 

h = j^~ E Mmj} (19) 

while the surface free energy density / s is defined as the excess from the average bulk 
free energy 

F = Y, M m i) = L Mb) = L n f b + 2f s . (20) 
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3. Magnetic properties 

3.1. Order parameter profile and critical temperature 

The order parameter profile is determined numerically by a Newton-Raphson method, 
starting with arbitrary values for the initial trial profile rrij. Equation (|]) provides a 
system of L coupled non-linear equations 



G i (m 1 ,m 2 , 



0. 



for the components of the vector m 
first order Taylor series: 



1,2, 



{m 1 ,m 2 , 



,L (21) 
mj, which can be expanded in a 



GAm + 5m) = GAm) + - — -5m j + 0(5m 
A set of linear equations follows for the corrections 5rh 

-GAm) 



9Gj 



(22) 



(23) 



which moves each function Gi closer to zero simultaneously This technique is known to 
provide a fast convergence towards the exact solution. Typical examples of the profile 
obtained for the Fibonacci perturbation are shown on figure |3|. 

The magnetization profile decreases as the temperature is increased and vanishes for 
some size- dependent effective value of the critical temperature (jl c (L) = 1 — (0 C (L)) _1 . 
This value may be obtained through a recursion relation deduced from the equation of 
state. In the high temperature phase, when h = 0, equation (Q) can be rewritten as a 
homogeneous system of linear equations: 



Grh 



V 



ai 


-1 








-1 


a 2 


-1 








-1 




-1 














( m 1 \ 

m 2 



(24) 



where atj 



2 + /i + rfj. If the determinant Dl(/j,) = Det G(/i) is not vanishing, the null 
vector m = provides the satisfying unique solution for the high temperature phase. 
The critical temperature is then defined by the limiting value fi c {L) which allows a non- 
vanishing solution for m, i.e. D^(yU c ) = 0. Because of the tridiagonal structure of the 
determinant, the following recursion relation holds, for any value of \x: 

D L (n) = a L D L _ x (y) - D L _ 2 (fi), (25a) 
DM = 1, (256) 
Dx(fi) = ai . (25c) 
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Figure 3. Order parameter profiles for a perturbation r = 2 and three values of the 
temperature below the critical point. The size of the chain is L = 144. 
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Thus we can obtain fi c {L) for different sizes L from 144 to 46368 and estimate the 
asymptotic critical point by an extrapolation to infinite size. This technique allows a 
determination of the critical temperature with an absolute accuracy in the range 10 -7 
to 10~ 9 depending on the value of the amplitude r. 



3.2. Surface and bulk spontaneous magnetization behaviours 

The boundary magnetization mi vanishes at the same temperature than the profile 
itself. First of all, the influence of finite size effects [^2| has to be studied. This is 
done by the determination of the profiles for different chains of lengths given by the 
successive sizes of the Fibonacci sequence L = 1,2, 3, 5, 8, 13, 21, 34 . . . The boundary 
and bulk magnetization in zero magnetic field are shown on figure ion a log-log scale. 




Figure 4. Log-Log plot of the bulk and boundary magnetization v.s. the reduced 
temperature t — fi c — fj, for two values of the aperiodic amplitude r and for different 
sizes of the chain from 144 to 46368. Finite-size effects occur when the curves deviate 
from the asymptotic straight line. The insert shows the behaviour of the magnetization 
with the temperature. 

The finite size effects appear in the deviation from the straight line asymptotic 
behaviour. These effects are not too sensitive, as it can be underscored by considering 
the deviation of the curve for a size L = 17711, which occurs around t = \i c — ji ~ 10 -7 , 
i.e. very close to the critical point. 
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The expected marginal behaviour is furthermore indicated by the variation of 
the slopes with the aperiodic modulation amplitude r and is more noticeable for the 
boundary magnetization than in the case of the bulk. 

A more detailed inspection of these curves also shows oscillations resulting from the 
discrete scale invariance of the system and the asymptotic magnetization can thus 
be written 

m(t) = t p fh{r v ) (26) 

where fh(t~ y ) is a log-periodic scaling function of its argument. We make use of this 
oscillating behaviour to obtain a more precise determination of the critical temperature 
(in the range 1CT 11 to 10~ 12 ) and of the values of the bulk and surface exponents by 
plotting the rescaled magnetization mt~^ as a function of \wt~ v as shown on figure [| in 
the case of the first layer. 




Figure 5. Periodic oscillations of the rescaled boundary magnetization m\t~^ v.s. 
lni~' y . The deviation from the oscillating behaviour for large values of the correlation 
length t~ u is due to finite-size effects. The insert shows the oscillations of the rescaled 
boundary magnetization for different values of r after substraction of a constant 
amplitude. 

The values of \i c and f3i that we consider suitable are the ones which allow an 
oscillating behaviour for the widest interval in the variable lnt~ u . A modification of 
the boundary exponent (3\ would change the average slope of the oscillating regime. 
This could be due to corrections to scaling, but, if such corrections really existed, they 
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should cancel in this range of temperatures (in the oscillating regime, t goes to values 
as small as 1CT 9 ). The other parameter, // c , modifies the number of oscillations and 
we have choosen a value leading to the largest number of such oscillations. A poor 
determination of the critical point fi' c = /x c + A/x c would indeed artificially introduce a 
correction to scaling term, since t@ = (f + A/i c )^ ~ t' 13 (1 + j3^r-). 



Table 2. Numerical values of the critical temperature and the magnetic exponents for 
the surface and bulk magnetizations. The figure in brackets gives the uncertainty on 





the last di 


=rit. 










r 






surface 




bulk 






0i 


0l 


0s 


0{n-l) 


0b 


.1 


.963977634341 (5) 


1.00036 (2) 




.0002 (2) 


.500087 (1) 


.5002 (2) 


.2 


.93187679929 (2) 


1.00146 (2) 


1.0015 (1) 




.50033 (1) 




.3 


.90314503363 (2) 


1.0034 (1) 


1.0034 (1) 




.50072 (6) 




.5 


.85404149087 (2) 


1.0092 (1) 




.0094 (2) 


.50187 (1) 


.505 (1) 


.8 


.796437160887 (5) 


1.02214 (2) 




.0216 (2) 


.50419 (2) 




1. 


.76600595095 (2) 


1.0327 (1) 




.0302 (2) 


.505777 (2) 


.516 (1) 


1.5 


.70902241601 (2) 


1.0621 (1) 






.50943 (1) 




2. 


.67010909237 (2) 


1.0913 (1) 




.087 (1) 


.51186 (3) 


.538 (1) 


2.5 


.64234629279 (2) 


1.1178 (1) 






.51294 (1) 




3. 


.621796760462 (5) 


1.1410 (1) 


1.1410 (1) 


.133 (1) 


.5132 (1) 


.555 (1) 


3.5 


.60610567508 (2) 


1.1602 (1) 






.51327 (4) 




4. 


.593804120472 (5) 


1.1766 (1) 




.1692 (4) 


.5130 (1) 


.563 (1) 


4.5 


.58394117369 (2) 


1.1904 (1) 






.5122 (1) 




5. 


.5758805295248 (5) 


1.2026 (1) 




.195 (1) 


.51125 (2) 


.567 (1) 



The corresponding values of C , Pi and /3( n -i) are given for several values of the 
perturbation amplitude r in table |[ The critical exponent associated to the right 
surface (mi) of the Fibonacci chain has also been computed for different values of r 
for the largest chain size. It gives, with a good accuracy, the same value as for the left 
surface (mi) as it can be seen by inspection in the table. The aperiodic sequence is 
indeed the same, seen from both ends, if we forget the last two digits. 

Furthermore, the profiles of figure ^ clearly show that the sites of the chain are 
not all equivalent and the magnetization profiles can be locally rescaled with different 
values of the exponents depending on the site |I8| . Thus, after the local quantities, the 
computation of the surface and mean bulk magnetizations enable us to determine the 
critical exponents respectively written (3 S and (5b and given in table |2|. 

From our values, one obviously recovers the usual unperturbed ordinary transition 
values of the exponents when the perturbation amplitude goes to zero. 
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3.3. Susceptibility and critical isotherm 

Taking account of a non vanishing bulk magnetic field in equations (Q) and (7), one 
can compute the magnetization in a finite field and then deduce the critical isotherms 
exponents 6r n -i) an d Si from the behaviours of the local magnetizations m( n _i) and mi 
with respect to h: 

m(n-i) ~ h 1/s ^\ rm ~ h l/Sl , t = 0. (27) 




Figure 6. Rescaled equations of state for the boundary and bulk magnetization 
for r = 2. The values of the temperature are 9 = 0.670090, 0.670094, 0.670097, 
0.670100, 0.670103, 0.670105 below C and 0.670111, 0.670113, 0.670115, 0.670118, 
0.670121, 0.670125 above 9 C - Top: scaling functions , the insert shows the boundary 
magnetization as a function of the bulk magnetic field. Bottom: same as above for the 
local bulk magnetization. 

This time, a direct log-log plot allows a precise determination of the exponents and 
the rescaled equation of state confirms the validity of the estimate since we obtain a 
good data collapse. In the case of the boundary magnetization, the scaling assumption 
takes the following form under rescaling by an arbitrary factor b: 

mi (t, h) = b- pi/v mi {b yt t, b yh h), (28) 
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where y t is given by the inverse of the correlation length exponent y t = 1/V and the 



value of the magnetic field anomalous dimension y^ follows the requirements of (27): 
Vh = PiSi/v = (38/ v. The choice b = t~ u for the rescaling factor then leads to a universal 
behaviour expressed in terms of a single scaled variable: 

m 1 (t,h) = t*f± L (hr A ) (29) 

where A = f3\8\ is the so-called gap exponent, is a universal scaling function and 
± refers to the two phases 9 > 9 C and 9 < 9 C . This may then be checked by a plot of 
m\t~^ x v.s. ht~ A shown on figure |6] and the same type of universal function have been 
obtained for the local bulk site m^.ijt - ^"- 1 ' = fm (n _ x) (ht~ A ). The values of 8i and 
<5( n _i) are given in table |3|. 



Table 3. Numerical values of the critical exponents associated to the critical isotherms 
and the susceptibilities. 7b and 5^, correspond to the behaviour of the mean bulk 
magnetization mj. The figure in brackets gives the uncertainty on the last digit. 

surface bulk 



r 


7i 


Si 




Is 




6s 


7(n-l) 


<5(n-l) 


76 


s b 


.1 


.5013 (2) 


1.5024 


(2) 


1.498 


(1) 




.9997 (1) 


2.9989 (1) 


1.0005 (1) 




.2 


.5006 (2) 


1.5004 


(2) 








.9993 (2) 


2.9972 (3) 






.3 


.4992 (2) 


1.4977 


(2) 








.9993 (2) 


2.9949 (9) 






.5 


.4958 (2) 


1.4901 


(2) 


1.493 


(1) 


312 (11) 


.9989 (2) 


2.9895 (9) 


.99790 (2) 


2.98136 (2) 


.8 


.487 (1) 


1.4751 


(3) 


1.486 


(1) 


85 (2) 


.9986 (3) 


2.981 (2) 






1. 


.4796 (2) 


1.4641 


(2) 


1.480 


(2) 


53 (1) 


.9985 (4) 


2.9744 (9) 


.99253 (2) 


2.93144 (2) 


1.5 


.4568 (2) 


1.4378 


(4) 








.9988 (7) 


2.963 (3) 






2. 


.4316 (2) 


1.4135 


(1) 


1.438 


(2) 


16.37 (2) 


.999 (2) 


2.9571 (1) 


.9792 (1) 


2.82375 (3) 


2.5 


.412 (1) 


1.3845 


(G) 








.9992 (9) 


2.954 (2) 






3. 


.388 (1) 


1.3484 


(G) 


1.394 


(2) 


11.2 (2) 


.9988 (6) 


2.952 (2) 


.9660 (1) 


2.7513 (2) 


3.5 


.372 (1) 


1.3108 


(5) 








.9986 (5) 


2.949 (2) 






4. 


.354 (1) 


1.2989 


(4) 


1.360 


(2) 


10.01 (5) 


.9988 (8) 


2.948 (2) 


.9619 (5) 


2.6976 (3) 


4.5 


.341 (1) 


1.2571 


(2) 








.999 (1) 


2.950 (2) 






5. 


.328 (1) 


1.2467 


(6) 


1.330 


(2) 


8.3 (4) 


.999 (2) 


2.953 (2) 


.9514 (2) 


2.65938 (2) 



The behaviours of m s and m;, with h at the critical point lead to the values of 5 S 
and 8b, also listed in table [| We can point out the low accuracy in the determination 
of 8 S since the slope of the log-log plot of m s v.s. h is quite small when r reaches the 
unperturbed value r = 0. 

The derivative of equation (|28|) with respect to the bulk magnetic field h defines the 
boundary susceptibility xi which diverges as the critical point is approached with an 
exponent 71 . Numerically, the boundary magnetization is calculated for several values of 
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the bulk magnetic field (of the order of 10~ 9 ), and \\ follows a finite difference derivation. 
The bulk local susceptibility X(n-i) ma ¥ be obtained in the same way. Log-periodic 
oscillations also occur in these quantities and the determination of the exponents can 
be done in the same way as in the previous section for the magnetization. Again, the 
accuracy of the result is confirmed by the rescaled curves for the susceptibilities, for 
example X(n-i)^ 7(n_1) = /x< n „ 1) (^~ A ) shown on figure ^exhibits a good data collapse 
on two universal curves for 9 < 9 r and 9 > 6 r . 




Figure 7. Rescaled bulk susceptibility giving the behaviour of the universal functions 
/± below and above 8 C for r = 2. The values of the temperature are the same as 
in figure 6. The inserts show the behaviours of X(n-i) &s a function of h for the same 
temperatures (left), and the singularities of both X(n-i) an d xi m zero magnetic field 
as a function of 9 (right). 

The values of the exponents are given in table |^ which presents also 7 S and 7&, 
associated to the surface and average bulk magnetization field derivatives. 

4. Specific heat 

According to the definitions of section [|, the surface and bulk free energies are also 
defined as follows: 

F S = ^(F FBC -F PBC ), (30 a) 

F b = F PBC . (306) 
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where F FBC and F PBC denote the total free energies of aperiodic chains with free and 
periodic boundary conditions respectively and are obtained numerically using equations 
(0) and ©. 

The expected singular behaviours of the free energy densities 

f s (t,h) = t 2 - a °f s (hr A ), (31a) 

f b (t,h) = t 2 -^f b (hr A ), (3ib) 

where the dependence of f s with the local magnetic surface field hi has been omitted 
since we always consider the case hi = 0, lead to the surface and bulk specific heat 
exponents. The values of a s and a b are simply deduced from the slopes of the log-log 
plots of f s and v.s. t. 

In figure we show the bulk free energy density amplitude fbt ab ~ 2 as a function 
of lnt _t/ for r = 2. It exhibits the same type of oscillating behaviour than the rescaled 
magnetisation of figure ||. 




Figure 8. Rescaled bulk free energy density fbt ab 2 v.s. bit v for r = 2. The 
amplitude of the bulk free energy density exhibits log-periodic oscillations. 



The surface and bulk specific heat exponents are collected in table f|. The bulk 
specific heat discontinuity of the homogeneous system is washed out in the perturbed 
system, since a b < 0. 



17 



Table 4. Numerical values of the specific heat critical exponents. The figure in 
brackets gives the uncertainty on the last digit. 



r 






.1 


0.5f496 (7) 


-0.00031 (1) 


.5 


0.50112 (5) 


-0.00733 (1) 


.8 


0.48448 (5) 


-0.01709 (1) 


1. 


0.47075 (4) 


-0.02462 (1) 


2. 


0.40265 (1) 


-0.05924 (1) 


3. 


0.35077 (4) 


-0.07813 (1) 


4. 


0.31516 (3) 


-0.08579 (1) 


5. 


0.28935 (6) 


-0.08805 (1) 



5. Discussion 



We have calculated numerically several surface and bulk critical exponents for a marginal 
aperiodic system within mean field theory. The marginal aperiodicity leads to exponents 
which vary continuously with the amplitude of the perturbation r. The variations of 
these exponents are shown on figure [| as a function of r. 

The comparison in table |2] between the bulk exponent /3& and the local one fl( n -i) 
clearly shows that it is no longer possible, in this aperiodic system, to define a unique 
bulk exponent, as it was already suggested by the possibility of a local rescaling of the 
profiles with position-dependent exponents which suggests a multiscaling behaviour. A 
constant value yt = 1/v is consistent with continuously varying exponents, in order 
to keep a vanishing crossover exponent which ensures that the marginality condition 
remains valid for any value of the aperiodicity amplitude r. For y^ on the other hand, 
there is no such reason. From this point of view, equations like ( f28|) are not exact 
since a unique field anomalous dimension y^ has no real significance. It follows that 
the universal functions in figure § and figure [7] only give an approximate picture of the 
scaling behaviour in this system, since they involve the gap exponent A = yh/yt- The 
good data collapse has to be credited to the weak variation of the exponents with the 
perturbation amplitude r. 

On the other hand, the scaling laws involving the dimension of the system are 
satisfied in mean field theory with a value of d equal to the upper critical dimension d*. 
As for the 2d Ising model with a marginal aperiodicity [[17], FJ], one expects a strongly 
anisotropic behaviour in the Gaussian model. It yields a continuous shift of the upper 
critical dimension with the pertubation amplitude, d*(r), since the value d* = 4 for a 
critical point in the homogeneous 4 theory follows Ginzburg's criterion for an isotropic 
behaviour. Hyperscaling relations should thus be satisfied for the mean field exponents 
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Figure 9. Variations of the surface and bulk exponents with the perturbation 

amplitude r (boundary exponents: - — , surface: — • — , local bulk: , mean 

bulk: - -). 



with d*(r): 



2 - a b = ud*(r), (32a) 
2-a s = u{d*(r)-l). (326) 

We can make use of these relations to obtain an estimate of the upper critical dimension 
d*(r) for this aperiodic system. The corresponding results are given in table g. 

The two determinations are in good agreement for small values of the perturbation 
amplitude. The discrepancy at larger values of r suggests that the precision in the 
determination of the exponents has probably been overestimated, but the variation of 
the upper critical dimension with the perturbation amplitude is clear and should be 
attributed to an anisotropic scaling behaviour in the corresponding Gaussian model. 

One can finally mention that a mean field approach for relevant aperiodic 
perturbations would be interesting. Many cases of aperiodic sequences with a wandering 
exponent u > — 1 are known, they constitute relevant perturbations in mean field theory. 
In the case of the 2d layered Ising model with relevant perturbations, a behaviour which 



looks like random systems behaviour, with essential singularities, was found flTB|], and 



19 



Table 5. Numerical values of the upper critical dimension d*(r) deduced from 
hyperscaling relations. 



r 


(2 - ab)/v 


(2 - a s )/v + 1 


.1 


4.00 


3.97 


.5 


4.01 


4.00 


.8 


4.03 


4.03 


1. 


4.05 


4.06 


2. 


4.12 


4.19 


3. 


4.16 


4.28 


4. 


4.17 


4.37 


5. 


4.18 


4.42 



the same type of situation can be expected within mean field approximation. 
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